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I. 



INTRODUCTION 



A body's resistance to flow is strongly affected by what 
surrounds it. When two bodies are in close proximity, not 
only the flow about the downstream body but also that about 
the upstream body may be influenced. Examples include con- 
denser tubes in heat transfer; columns in pressure suppression 
pools of nuclear reactors; risers, piles, and other tubular 
structures in off-shore engineering; turbine and compressor 
blades in mechanical or aerospace engineering; and high-rise 
buildings, cooling towers, and transmission lines in civil 
engineering. The quantification of the interference effects 
in terms of the pressure distribution, lift and drag forces 
on individual members, vortex shedding frequency, and the 
dynamic response of members of the array in terms of the 
governing flow and structural parameters constitute the 
essence of the problem. 

There are infinite numbers of possible arrangements of 
two or more bodies positioned at right or oblique angles to 
the approaching flow direction. In wavy or time-dependent 
flows one needs the lift, drag and inertia coefficients for 
all members of the array. Evidently, the members of the 
array may not be all parallel and normal to the flow. The 
quantification of the flow interference on lift, drag and 
inertia coefficients for cylinders with relative inclinations 
and spacings in a design-wave environment is an exceedingly 
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complex problem. In the absence of data on the interaction 
between drag and inertia coefficients in wavy and harmonic 
flows for cylinder groups one is tempted to use steady-flow 
results for the drag coefficient and the unseparated potential 
flow results for the inertia coefficient. The inertia coef- 
ficient for a group of cylinders in inviscid unseparated 
flow may be obtained through the use of the method of images 
/~ 1, 2_7 or through the use of the linear potential theory 
including wave diffraction / — 3 , 4_/. These analyses do not 
deal with the effects of separation and vortex shedding. Con- 
sequently, the results are more appropriate to the determina- 
tion of earthquake forces and wave forces on large bodies 
rather than to the evaluation of the inertial components of 
the force in the drag/inertia dominated regime. 

It is evident from the foregoing that experiments must be 
carried out with carefully selected tube bundles and cylin- 
der combinations in order to develop some understanding of 
the flow interference and to provide data for body combina- 
tions of practical importance. 
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II. REVIEW OF PREVIOUS INVESTIGATIONS 



A careful review of flow interference between two circular 
cylinders in various arrangements in steady flow has been pre- 
sented by Zdravkovich / 5_7 where an extensive list of re- 
ferences may be found. Numerous studies have shown that the 
changes in drag, lift and vortex shedding are not necessarily 
continuous. In fact, the occurrence of a fairly abrupt change 
in one or all flow characteristics at a critical spacing is 
one of the fundamental observations of flow interference in 
cylinder arrays. 

For the tandem arrangement (one cylinder behind the other, 
see Figure 1) , it has been shown that at relative spacings 
S/D<3.5 there is a strong mutual interference between the two 
cylinders. This critical distance increases with the bluff ness 
of the bodies. For two plates in tandem, the critical spacings 
is about 10 times the plate width / 6 /. In general, the 
tandem arrangement has a strong effect on the drag and is sen- 
sitive to spacing. The upstream cylinder takes the brunt of 
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Fig. 1 Tandem arrangement of two cylinders 
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the burden and the total drag for the group is smaller than 
the sum of the drag forces acting on each cylinder in isola- 
tion in a tunnel with the same blockage. The drag coefficient 
for the downstream cylinders is reduced partly by shielding 
and partly by the occurrence of earlier transition in the 
boundary layers due to the turbulence in the approaching flow. 



FLOW 



-H D |^- 




Fig. 2 Side-by-side arrangement of two cylinders 



The side-by-side arrangement (see Figure 2) exhibits a 
discontinuity in the flow and measured forces for spacing 
ratios smaller than a critical value (S/D<2.2 for two circu- 
lar cylinders and about 4 two plates) . The observed dis- 
continuity (a switching phenomenon and bistable nature of 
the flow between plates) results in two values of the drag 
coefficient. This results from the mutual interference of 
the vortices on the adjacent sides of the vortex streets. 

The bistable nature of the flow may be avoided by increasing 
the size of one of the cylinders. The sum of the bistable 
high and low drag is often less than twice the drag of the 
single cylinder. For two plates in side-by-side arrangement. 
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the said sum may be about 10 percent larger than twice the 
drag of the single plate. This increase strongly depends on 
the width of the flow field relative to the total width of 
the bodies, i.e., on the blockage ratio. For relatively 
large blockage, the flow is forced through the openings 
between the cylinders rather than around the total configura- 
tion. Consequently, the use of ordinary shielding and block- 
age factors for groups of cylinders is meaningless. Experi- 
ments must be conducted in channels or tunnels with very small 
blockage in order to obtain valid force-transfer coefficients. 
There is at present no means to separate the flow interference 
effects from the blockage effects for groups of cylinders. 

Additional work with tube arrays in steady flows has been 
reported by a number of people / 7 - 14 /. Ross / 15 7 con- 
ducted large scale wave tank tests for the case of one cylin- 
der on each side of the test cylinder in the range of critical 
Reynolds numbers. His results appear to indicate that the 
wave force increases significantly only when the spacing 
between two cylinders is less than about one diameter. 

The foregoing studies, conducted mostly in steady flows, 
show that the sum of the individual drag forces for isolated 
cylinders is often greater than the total force acting on the 
array as long as the local boundaries do not constrain the 
flow. They further show that the results obtained from the 
tandem and side-by-side arrangements cannot be generalized to 
predict the combined interference resulting from the staggered 
arrangement. 
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Gibson and Wang / 16_/ carried out two different experi- 
ments to determine the added mass of a series of tube bundles. 
The bundles consisted of tubes of uniform diameter d arranged 
either in a square configuration or a circular configuration. 

In the first series of experiments, they towed the model of a 
pile cluster under linear acceleration. In the second series, 
they have vibrated the model at its own natural frequency. For 
both cases , they have calculated the added mass through the 
use of the measured force and acceleration and plotted them as 
a function of the 'solidification ratio' defined by Zd^/irD 
where D is the pitch diameter of the bundle. Their results 
have shown that the added mass increases sharply after the 
solidification ratio reaches the value of 0.4 to 0.5. Beyond 
this value, the volume enclosed rather than the volume dis- 
placed by the structure becomes important. This result is 
disputed / 17_/ and the results of both series of tests are 
no more applicable to separated wavy or oscillatary flows 
about tube bundles than those predicted from the potential 
theory with or without diffraction effects. 

Relatively few studies have been carried out with oscil- 
lating tube bundles / — 18, 19_7* Bushnell / — 19_/ determined 
the interference effects on the drag and transverse forces 
acting on a single member of a two-cylinder configuration 
through the use of a pulsating water tunnel. He did not 
evaluate the drag and inertia coefficients through the use 
of a suitable method, e.g. Fourier averaging. Instead, he 
picked out the maximum force values which occurred in each 
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half cycle and averaged them over ten consecutive values so as 
to obtain a mean maximum force for each flow direction. The 
results have shown that the presence of neighboring cylinders 
significantly affects the forces on an individual cylinder of 
an array and the interference effect increases with increasing 
relative flow displacement. The maximum drag force on shielded 
cylinders was reduced relative to an exposed cylinder by up to 
50 percent. Bushnell has suggested that a design using a high 
Reynolds number single-cylinder drag coefficient applied through- 
out the array would have an extra margin of safety against maxi- 
mum drag loading due to interference effects. The transverse 
force could be 3 to 4 times larger for interior array positions 
than that for a single cylinder. Thus, a cylinder array, such 
as a riser, supported at regular intervals may exhibit very 
complex dynamic behavior. Some of the members between supports 
may undergo violent transverse oscillations. 

It is evident from the foregoing that there is very little 
systematic data for tube bundles and two cylinder interference. 

In view of this fact is was decided to carry out experiments 
with two bundle configuration consisting of 12 and 9 outer 
tubes, respectively. Subsequent to this investigation a 
detailed study of the interference effects of two cylinders 
in harmonic flow was undertaken. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 



A. OSCILLATING FLOW TUNNEL 

Experiments were carried out in a large U-shaped oscillating 
water tunnel (see Figure 3) . It has been used extensively at 
this institution over the past 6 years /” 20 Only salient 
features, most recent modifications, as well as the adaptation 
for this work, are briefly described here. 

The length of the tunnel has been increased from 30 feet to 
35 feet and its height from 16 feet to 24 feet. Previously, a 
butterfly-valve arrangement at the top of one of the legs was 
used to initiate the oscillations. In 1978, the tunnel was 
modified so that the oscillations could be generated and 
maintained indefinitely at the desired amplitude. For this 
purpose the output of a 2 HP fan was connected to the top of 
one of the legs of the tunnel with a large pipe. A small 
butterfly valve, placed in a special housing between the top 
of the tunnel and the supply line, oscillated harmonically 
at a frequency equal to the natural frequency of the oscilla- 
tions in the tunnel. The oscillation of the valve was per- 
fectly synchronized with that of the flow through the use of 
a feedback control system. The output of a pressure trans- 
ducer (sensing the instantaneous acceleration of the flow) 
was connected to an electronic speed-control unit coupled to 
a DC motor oscillating the valve plate. The circuit main- 
tained the period of oscillations of the valve within 0.001 
seconds. The fluid oscillated with a period of 6.000 seconds. 
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Fig. 3 U-shaped oscillating water tunnel 



The amplitude of oscillations was varied by constricting or 
enlarging an orifice at the exit of the fan. The flow oscil- 
lated at a given amplitude as long as desired. 

Two identical force transducers, one at each end of the 
cylinder, were used to measure the instantaneous in-line and 
transverse forces. A special housing was built for each gage 
so that it can be mounted on the tunnel window and rotated to 
measure either the in-line or the transverse force alone. 

The calibration of the each gage was accomplished by 
hanging loads in the middle of the cylinder after setting both 
gages to sense only the transverse force (here in the vertical 
direction) . Repeated calibrations have shown that (a) the 
gages were absolutely linear up to 500 pounds; (b) the gages 
yielded the same signal for loads applied either upward or 
downward; and that (c) the gages, together with the electronic 
system to which they are attached, were capable of sensing 
loads as small as 0.02 pounds. 

The displacement, velocity, and acceleration of the fluid 
in the tunnel were obtained through the use of several methods 
described in / 20 /• 

The data were recorded in analog form and then digitized 
and processed through the use of an HP9845 computer-digitizer 
system. 

B. TEST BODIES 

The tube bundles were made of aluminum rods and mounted 
between two circular stainless steel plates . Two configura- 
tions were chosen for test purposes. These configurations 
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Fig. 4b Riser 2 configuration 
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roughly correspond to the most commonly used riser shapes in 
the off-shore industry. Figures 4a and 4b show the two riser 
configurations together with their characteristic dimensions. 

The two cylinder configuration was arranged in a similar 
manner (see Figure 5) . The active cylinder was mounted on the 



force transducer as described previously. The passive cylin- 
der spanned the entire test section and was rigidly mounted 
on the tunnel walls. Both cylinders were smooth and 3 inches 
in diameter. The angle a was varied from zero to 90 degrees 
in 30 degree increments. The N in cylinder spacing ND was 
varied from 1.5 to 3.5. The in-line and transverse forces 
acting on the active cylinder were measured and analyzed in 
terms of the appropriate lift coefficients and the drag and 
inertia coefficients through the use of the Fourier analysis. 
The details of the definition of the various force transfer 
coefficients are described in the following. 




Fig. 5 General arrangement of two cylinders 
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IV. FORCE TRANSFER COEFFICIENTS 



A. TUBE BUNDLES 

The in-line force, based on Morison's equation is written 
as 



2 i i 

F = -0.5/C,U LED. Cosut Cosut 
dm 1 1 1 

( 1 ) 

+ 0.25tt/LC E D?U cuSinout 
m l m 



in which U^ is the maximum velocity in a given cycle. 
Denoting 



D = E D. 
a l 




and inserting in equation (1) , one has 



F 

0.5/D LU 2 
e m 



D 

-C^* (- — ) • | Cos wt | *Cos cut 
e 



+ 



7 T 



U T 
(t?-) 



C ’ Sincut 
m 



( 2 ) 



in which T = 2tt/u . 

For an oscillating flow represented by U = -U^ coscut the 
averages of C^ and C^ are given by / 20_/ 
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T D ~ F Cosojt 

c. = | • <=A / 2l1 -S _ 

d 4 D a 0 /LD U 2 

e m 



d at 



(3) 



and 



2U T 

C = / 

TT D O 



~ F *Sino)t 
2 tt m 

? L D U 2 
e m 



d at 



(4) 



in which represents the measured force. 
The lift coefficient is defined as 



C L " 



0.5/LD U 
e m 



(5) 



in which F L represents the amplitude of the transverse force 
and L, the length of the cylinder. 

The vortex shedding frequency for the entire bundle was 
expressed in terms of the Strouhal number defined by 



St = 



f D f T 
v e v 



U 



m 



K 



in which f represents the frequency of vortex shedding; T, 
the period of flow oscillation; and K, the Keulegan-Carpenter 
number defined by 



K = 



U T 
m 

D 
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B. TOO CYLINDER ARRANGEMENT 

The in-line force acting on the active cylinder is not in- 
dependent of the direction of the fluid motion because of the 
interference effect of the passive cylinder. In the interval 
ir/2 <ojt< 3 it/ 2 the direction of flow is opposite to that in the 
interval 3 tt/ 2 <ajt <5ir/2 . With this fact in mind Morison's equa- 
tion was written as 

2 

j 7 =- • C sinurt - C, | Cos cot | Co suit (6) 

0.5/D U m m 
m 



and C and C, were determined from 
m d 



4 -U -T , /0 _ . 

m , 3 it/ 2 F • Smoot , 

mi tt 3 F tt/2 /U 2 LD 

m 



( at) 



„ _ 3 r 3 tt/2 F-Coscut A t.*.\ 

C,, =-? / — ^ * d(wt) 

d4 Z tt/2 / U Z LD 

m 



'm2 



4U T 

m r 5ir/2 F.Sincot ,, .. 
“XT / 3n/2 ZZ1 ' d(Ut) 



TT D 



AJ LD 
m 



and 



M2 



3 r 5ir/2 F-Coscot .. 
- j f 5 d(cot) 

3 tt/2 /U LD 
m 



(7) 



(8) 



(9) 



(10) 
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in which C 



and represent the inertia and drag coefficients 

for the interval tt/2 curt <3 tt/2 ; and C and / the corresponding 
coefficients for the interval 3tt/2 <ajt <5tt/2 . In general these 
two sets of in-line force coefficients are not identical accept 
when- a=90°. 

The evaluation of the transverse-force coefficients re- 
quires a distinction not only in the direction of flow but 
also in the direction of the transverse force. As the flow 
proceeds in a given direction the lift force acting on the 
active cylinder does not remain symmetrical primarily because 
of the proximity of the passive cylinder, situated asymmetri- 
cally relative to the direction of the flow. It is also re- 
cognized that the random nature of the vortex shedding does 
not allow a perfectly symmetric lift force even in the absence 
of any interference effects. Consequently, one must define 
four lift coefficients. In the interval tt/2 < at <3 tt/2 , one 
has (see Fig. 6a) 



Ll 



'Ll 



0.5/U DL 
m 



C 



L2 



L2 



0 .5/-IT 



m 



D-L 



and in the interval 3tt/2 <cot <5ir/2 , (see Fig. 6b) 
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and 



where 



L3 

'L3 ~ 2 

LJ 0.5/r DL 
M 



L4 



L4 0.5/U 2 DL 
m 



F ^ through ^ are shown in Figures 6a and 



6b. 
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Fig. 6a Definition sketch for lift forces 
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Fig. 6b Definition sketch for lift forces 
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The Strouhal number was defined in a manner similar to that 



previously described and evaluated through the use of 



St. = 
1 



f . -D 
vi 

U 

m 



in which i = 1, .. .4, indicating the lift force directions 

shown in Figures 6a and 6b. 

C. GOVERNING PARAMETERS 

It has been shown previously / 20_7 that the Fourier- 
averaged lift, drag, inertia coefficients for a smooth cylin- 
der immersed in a harmonically oscillating flow depend on 
the Reynolds number and the Keulegan-Carpenter number. In 
other words 



(C . , C , C } 
d m L 



= F. 



l 



U D U T 
m m . 

v 7 D 1 



In terms of a frequency parameter 6 defined by 



U D U T 

6 = (-E-) / (-§-) 




VT 
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one has 



{C 



d 



C _ / C x > 

m L 




U T 
m 

D 



) 



For the two cylinder arrangement one has to introduce two 
additional parameters namely, the angle a of the line joining 
the center of the two cylinders and the relative distance 
between the cylinders. Thus, all the force-transfer coeffi- 
cients cited previously depend on the Keulegan-Carpenter 
number, Reynolds number (or the frequency parameter) , the 
angle a, and the relative distance factor N (see Figure 5) . 
Thus one has 



C^(a coefficient) = F^(K,g,a,N) 
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V. DISCUSSION OF RESULTS 



A. TUBE BUNDLES 

The drag and inertia coefficients for the two riser con- 
figurations are shown in Figures 7a and 7b as a function of 
the Keulegan-Carpenter number for constant values of 6. The 
drag coefficient decreases gradually with increasing K. The 
inertia coefficient increases with increasing K and reaches 
values as large as 6. The use of the potential flow theory 
for unseparated flow yielded / 21 / an average inertia 

coefficient C* =2.17 for the two configurations. C* was 
m m 

defined by 



C* = (EC .D 2 )/ZD 2 
m mi i i 



The comparison of C*^ with that obtained experimentally has 
shown that some fluid mass is entrapped within the array and 
that neither the potential theory nor the diffraction analy- 
sis can adequately describe the behavior of the complex 
separated flow through the bundle. 

A comparison of the Figures 7a and 7b shows that the drag 
and inertia coefficents as defined here are nearly indepen- 
dent of the particular arrangement and the number of tubes 
comprising the bundle. 
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Experiments with different tube arrangements / — 21_7 have 
shown that the force coefficients are independent of the 
Reynolds number within the range of Reynolds and Keulegan- 
Carpenter numbers reported herein. 

The reason for the dependence of and on K is thought 
to be the dependence of the interaction of the wakes of the 
outer and inner pipes. The vortices in the wake of a given 
cylinder loose about 70 percent of their strength within 10 
cylinder diameters. Thus, for small values of K, the vortices 
generated by a small tube at the center front of the bundle 
arrive at the central tube as weak vortices. Consequently, 
each tube behaves more or less as if it were independently 
subjected to a turbulent harmonic flow. As K increases, not 
only the turbulence level but also the interaction between 
the wakes of the various cylinders increases. There is a 
certain amplitude of oscillation beyond which neither the 
interaction of the wakes nor the increase of the turbulence 
level affects the overall force acting on the bundle. A com- 
parison of the total drag force acting on the bundle with the 
sum of the drag forces acting on each cylinder in isolation 
in harmonic flow (at the corresponding K and Re values 
appropriate to each tube) shows that the former is about 10 
percent smaller. 

The lift coefficients, as given by Eqn. (5), for the two 
riser configurations are shown in Figures 7c and 7d. Apparently, 
both configurations yield nearly identical lift coefficients. 

As will be seen later the rapid decrease of the lift coefficient 
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with increasing K is in conformity with that for a single 
cylinder. As the Keulegan-Carpenter number increases, each 
member of the bundle sheds an increasing number of vortices 
with different strengths and frequencies. 

Furthermore the shielding or proximity effect of the tubes 
results in a highly confused vortex system. This in turn 
decreases the coherence of the regular vortex shedding, re- 
sulting in very low lift coefficients. Evidently, a greater 
understanding of the lift force requires the determination of 
the lift force acting on each member individually and the 
consequence of tube proximity. 



B. DRAG AND INERTIA COEFFICIENTS FOR THE TWO CYLINDER 
ARRANGEMENT 

The drag and inertia coefficients for N=1.5, N=2, N=2.5, 

and N=3.5 (only for a=0) and for a=0°, a=30°, a=60°, and a=90° 

are shown in Figures 8 through 20. Each figure represents 

either the drag or the inertia coefficient for two different 

flow directions. For example, Figure 8a shows the drag 

coefficients C,.. and C,~ for a=0° and N=1.5. The coeffi- 
dl d2 

cients C^ is for the flow directed from passive towards 
active cylinder. The coefficient C^ is f° r the flow dir- 
ected from the active towards the passive cylinder. In 
addition, included in each figure is the drag or inertia 
coefficient for the single cylinder of the same size. 

Figure 8b shows the inertia coefficients C ^ and C m9 for 



'm2 



a=0 and N=1.5. The coefficients C . and C are related to 

ml m2 

the flow direction in the same manner as the coefficients 



C dl and C d2 * 
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Figures 8a, 9a, 10a, and 11a show that for the tandem 
arrangement corresponding to a=0° the drag coefficient C ^ 
is always larger than . The difference between the two 
drag coefficients decreases with increasing N (larger spacing) . 
This finding is in conformity with the measurements reported 
by Zdravkovich / 5_/ for the tandem arrangement in steady 
flow. For relatively small spacings (smaller N) the wake of 
the passive cylinder is highly turbulent. This causes earlier 
separation on the active cylinder and hence in smaller 
values . One may also express the same fact by stating that 
the active cylinder is shielded by the passive cylinder when 
the flow sees the passive cylinder first. As the spacing 
between the two cylinders increases, and C ^ gradually 

approach each other and single cylinder value . Even for 
N=3.5, C ^2 is slightly larger than because of the afore- 

mentioned effect of turbulence in the wake of passive cylinder. 

The inertia coefficients C ^ and for the tandem arrange- 

ment are shown in Figures 8b, 9b, 10b, and lib. Unlike the 
drag coefficient, the inertia coefficient C £ is always 
smaller than C ^ . This is not entirely unexpected. All drag 
and inertia evaluations for single cylinders in waves and 
harmonic flows have shown that / 20_/ the drag and inertia 
coefficients follow opposite trends. In other words, when 
the drag coefficient increases the inertia coefficient de- 
creases and vice versa. The reason for this is embedded in 
the formulation of the Morison's equation (See Eq. (1)). When 
the phase angle between the velocity and the maximum force 
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is small the total force is primarily due to the velocity- 
square dependent drag and the contribution of the inertial 
force is negligible. In fact such flows are termed drag- 
dominated flows. When the phase angle between the velocity 
and maximum force is large, the total force is primarily due 
to the acceleration-dependent inertia. Such flows are termed 
inertia-dominated flows . Thus , it is apparent that the drag 
and inertia coefficients should follow opposite trends. 

As the relative spacing increases (See Fig. lib) C ^ and 
C m 2 approach each other. However, the difference between two 
coefficients remains relatively larger than that between 
and C^ 2 • The experimental fact is that the spacing between 
the two cylinders in tandem arrangements must be considerably 
larger than 3.5 diameters for the inertia coefficients to 
become independent of the proximity effects. It is also clear 
that structures in close proximity can cause relatively more 
significant changes in the inertia coefficient than in the 
drag coefficient. 

Figures 8b, 9b, 10b, and lib also show that the inertia 
coefficient for single cylinder (dashed lines) is nearly the 
same as C 2 • In other words , when the flow sees the active 
cylinder first the proximity effect on the inertia coefficient 
C m 2 is negligible. It is also apparent that when the flow 
sees the passive cylinder first the inertia coefficient C ^ 
is considerably increased. The important conclusion to be 
reached on the basis of these results is that the inertia 
coefficient of a shielded cylinder may increase considerably 
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and the drag coefficient may decrease significantly. 

The other extreme of the two cylinder arrangement is the 
side-by-side arrangement. The discussion of this case will 
be taken up next since the force coefficients for all other 
intermediate values of a follow trends intermediate to the 
tandem and side-by-side arrangements. 

Figures 18a, 19a, and 20a show and for N=1.5, N=2 , 

N=2 . 5 for the side-by-side arrangement. Several facts are 
immediately apparent from these figures. For spacings as small 
as 1.5 diameter the drag coefficients are not greatly in- 
cluenced by the proximity effects. Secondly, and C ^ 

are very close to that for the single cylinder. The slight 
difference in the lower range of Keulegan-Carpenter numbers 
may be attributed to the increased blockage in the tunnel and 
experimental uncertainty which was found to be about 5 per- 
cent. Experiments / — 5 _7 with steady flows about the same 
type of two-cylinder arrangement have also shown that the 
proximity effects are negligible for spacings larger than 
about 1.5 diameter. 

The inertia coefficients C ^ and C £ shown in Figures 18b, 
19b, 20b exhibit a somewhat different behavior. Even though 
they are nearly identical for any given spacing they are con- 
siderably larger than that for the single cylinder. Only for 
N=2.5 that they approach the single cylinder values. 

The two other arrangements intermediate to those discussed 
above behave as would be expected. Both the drag and inertia 
coefficients shown in Figures 12 through 17 vary between the 
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cases for a=0° and a=90° and both gradually approach the 
corresponding single cylinder values. The conclusions re- 
garding the effect of shielding remain valid for the inter- 
mediate arrangements. 



C. LIFT COEFFICIENTS FOR THE TWO CYLINDER ARRANGEMENT 

As in the case of drag and inertia coefficients, first the 
tandem arrangement, then the side-by-side arrangement, and 
finally, the intermediate arrangements will be discussed. 

For the tandem arrangement C L ]_ =C L 2 anc ^ C l 3 =C L4 ^ See F: “-9 S * 

6a and 6b) . The coefficients C^ or C^ are shown in Figures 
21-24 and the coefficients C^ or C T ^ are shown in Figures 
25-28. Also shown in each figure is the lift coefficient for 
a single cylinder of the same size. 

The data showed that the lift coefficient is strongly 
affected by the cylinder proximity, irrespective of the 
direction of the flow. For Keulegan-Carpenter numbers larger 
than about 20 the lift coefficients C ^ and C ^ remain prac- 
tically constant. Only for very small Keulegan-Carpenter 
numbers that C^ and increase rapidly, in conformity with 

that for a single cylinder. The reason for this is that for 
Keulegan-Carpenter numbers near 10, at most two vortices are 
shed during a given half-cycle. These vortices do not travel 
1.6 diameters in one direction as seen from the simple relation 
that K=2 tta/D where A is the amplitude of the flow oscillation. 
Consequently, the vortices shed from the active cylinder are 
only partially blocked by the passive cylinder for K smaller 
than about 10. Nevertheless, the effect of cylinder proximity 
is ever present as seen from a comparison of the data with 
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that for the single cylinder. As the spacing increases, the 
lift coefficients gradually approach that of the single cylin- 
der as seen in the Figures 24 and 28. 

In the case of side-by-side arrangement C =C , and 

LI Lj 

C L 2 =c L 4 * The coefficients or are shown in Figures 29- 

31 and the coefficients C ^ or C L 4 ' in Figures 32-34. The 
lift coefficients for this particular arrangement have been 
measured for only N=1.5, N=2, and N=2 . 5 . 

As the flow oscillates about the cylinders arranged in 
side-by-side position, the blockage introduced by the proximity 
of the cylinders forces the fluid to flow above the top and 
below the bottom cylinder. In other words, the velocity on 
the extreme sides of the passive and active cylinder increases 
while the velocity between the two cylinders decreases. This 
in turn increases the pressure on the sides of the cylinders 
nearest each other and decreases the pressure on the opposite 
sides. Consequently, or should be smaller than 

or C^. A comparison of the Figures 29-31 with Figures 
32-34 shows that C ^ or C L3 is in fact, always smaller than 
C L 2 or c L 4 • As the relative spacing increases all four 
coefficients gradually approach that for a single cylinder as 
seen from Figures 31 and 34. These results are in conformity 
with those obtained with steady flow about a side-by-side 
arrangement of circular cylinders. / 5 / 

For arbitrary arrangements one must deal with four dif- 
ferent lift coefficients, namely, C^/ ^ , and C T 



'L4 



The coefficient is shown in Figures 35-37 for a=30 and 
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for N=1.5, N=2 , and N=2 . 5 . Similarly; C 2 is shown Figures 
38-40; in Figures 41-43; and finally, in Figures 

44-46. The corresponding lift coefficients for a=60° arrange- 
ment is shown in Figures 47-58. 

A comparison of Figures 35-37 with Figures 38-40 shows 
that C ^ is always larger than • The reason for this is 
that the vortices shed from the passive cylinder pass over 
the top of the active cylinder and thereby reduce the pressure 
at the top of the active cylinder. The coefficients C ^ and 
C L 2 gradually approach that for a single cylinder as the 
spacing increases. 

When the flow first sees the active cylinder, C ^ remains 
relatively larger than C ^ because of the presence of the 
vortices near the top of the active cylinder which were shed 
during the previous half cycle. As the spacing between the 
two cylinders increases both and C l4 rapidly approach the 

lift coefficient for a single cylinder. 

For a=60° arrangement all four lift coefficients are 
practically unaffected by the cylinder proximity even for a 
spacing as small as 1.5 diameters and they are nearly equal 
to that of an isolated cylinder. Evidently, the effect of 
cylinder-cylinder proximity on the lift coefficient for a 
harmonically oscillating flow is confined to relatively small 
spacings and angular positions. 

The frequency of the lift force oscillations has been 
expressed in terms of the Strouhal number and evaluated for 
all cylinder arrangements. The results have shown that the 
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Strouhal number remains constant at an average value of about 
0.2. The only exception to that was the case corresponding 
to a=30° and N=1.5. 

Figures 59 and 60 show the Strouhal number as a function 
of K for a=0° and N=1.5. As N increased the Strouhal numbers 
for the two flow directions rapidly approached each other and 
a mean value of 0.2. 

Figures 61 and 62 show the Strouhal number for a=90° and 
N=2. Once again the lift coefficients for the two directions 
of flow are identical and gradually approach a value of about 
0.2. In the case of a=30°, the wakes interfere in an asymme- 
tric manner and the Strouhal number shows a strong dependence 
on the Keulegan-Carpenter number. Figures 63-68 show the 
Strouhal number as a function of K for various values of N 
for a=30°. Apparently, the Strouhal number decreases from 
values as large as 0.4 to about 0.15 as K increases. For 
spacings larger than about two diameters Strouhal number 
decreases to about 0.2. 
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VI. CONCLUSIONS 



An extensive investigation of the forces acting on tube 
bundles and arrangements of two cylinders has warranted the 
following conclusions: 

1. Experiments with two particular riser configurations 
have shown that the inertia coefficient is considerably 
larger than that predicted by the potential theory 

and shows that some fluid mass is entrapped within the 
bundle as a consequence of the solidification of the 
tube configuration. 

2 . For tandem arrangement of two cylinders the drag and 
inertia coefficients depend on both the Keulegan- 
Carpenter number and the relative spacing of the 
cylinders. As the amplitude of flow oscillation be- 
comes comparable or smaller than the gap between the 
two cylinders, the drag and inertia coefficients 
gradually approach those corresponding to an isolated 
cylinder, i.e., the wakes of the two cylinders do not 
interfere with each other. 

3. The drag and inertia coefficients for the side-by-side 
arrangement exhibit a similar behavior. For a relative 
spacing larger than about 2.5 diameters, the cylinders 
behave as if they were independent. 

4. The quantification of the lift force for an arbitrary 
arrangement requires four different lift coefficients. 
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depending on the direction of the flow and the direc- 
tion of the force. These coefficients have been eval- 
uated and it has been shown that the effect of cylinder- 
cylinder proximity is confined to relatively small 
spacings and angular positions. 

5. The Strouhal number for practically all two-cylinder 
arrangements remained nearly constant at about 0.2. 

Only in the case of 30° orientation did the strong 
asymmetric interference of the cylinder wakes result 
in a Keulegan-Carpenter-number dependent Strouhal 
number, ranging from 0.4 to 0.2. However, as the 
spacing increased the Strouhal number again approached 
a mean value of 0.2. 
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Fig. 7a C and C. versus K for riser 
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Fig. 7b C and C, versus K for riser 
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Fig. 7c C vs K for riser 



CD 



** 



CD 



CU 



n 000 1 



001 



eg 



01 



s> 



44 



Fig. 7d C vs K for riser 
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Fig. 8a C , and C, ~ versus K for N=1.5 and a=0° 
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Fig. 8b C , and C _ versus K for N-1.5 and 

ml mO 
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Fig. 9a C, and r,- versus K for N=2 and 
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Fig. 9b C , and C „ versus K for N=2 and 
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Fig. 10a C,. and C, 9 versus K for N=2.5 and 
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Fig. 10b C and C „ versus K for N=2.5 and 
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Fig. 11a C , , and C , 0 versus K for N=3 . 5 and 
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Fig. lib C . and C „ versus K for N=3.5 and 
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Fig. 12a C , and C,o versus K for N=1 . 5 and a=30 
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Fig. 12b C , and C 0 versus K for N=1.5 and c*=30 
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Fig. 13a C, and C, 0 versus K for N=2 and 
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Fig. 13b C and C „ versus K for N=2 and a=30 
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Fig. 14a C,, and C, 9 versus K for N=2.5 and a=30 
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Fig. 14b C , and C 0 versus K for N=2 . 5 and a=30 
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Fig. 15a C. and C,„ versus K for N=1 . 5 and a=60 
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Fig. 15b C and C „ versus K for N=1 . 5 and a=60 
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Fig. 16a C. 1 and C, 9 versus K for N=2 and a=60 
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Fig. 16b C , and C „ versus K for N=2 and a=60 
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Fig. 17a C, and C, versus K for N=2.5 and a=60 
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Fig. 17b C . and C 0 versus K for N-2.5 and a=60 
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Fig. 18a C, , and C „ versus K for N=1.5 and a=90 
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Fig. 18b C-,, and C-, 0 versus K for N=1.5 and a=90 
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Fig. 19a C,. and C, 0 versus K for N=2 and a=90 
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Fig. 19b C , and C 9 versus K for N=2 and a=90 
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Fig. 20a C,, and C, „ versus K for N=2.5 and a=90 
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Fig. 20b C and C ' versus K for N=2.5 and a=90 
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Fig. 21 C , or C 0 vs K for N=1.5 and 
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Fig. 2 2 C or C , vs I< for N=2 and 
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or C_ „ vs K for N=2 . 5 and 
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Fig. 24 C . or C T „ vs K for N=3.5 and 



Single Cylinder 
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Fig. 25 C _ or C . vs K for N=1 . 5 and 
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Fig. 26 C . or C vs K for N=2 and 
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Fig. 27 C or C . vs K for N=2 . 5 and 



Single Cylinder 
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Fig. 28 C - or C vs K for N=3.5 and 



Single Cylinder 
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Fig. 29 C, . or C vs K for N=1.5 and a=90 
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Fig. 30 C T , or C. vs K for N=2 and a=90 
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Fig. 31 C , or C vs K for N=2 . 5 and a=90 
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Fig. 32 C ~ or C vs K for N=1.5 and a=90 
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Fig. 33 C T 9 or C . vs K for N=2 and a=90 
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Fig. 34 C T ~ or C, . vs K for N=2.5 and a=90 
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Fig. 35 C T , vs K for N=1.5 and a=30 
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Fig. 36 C . vs K for N=2 a=30 
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Fig. 37 C . vs K for N=2.5 and 



Single Cylinder 
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Fig. 37 C - vs K for N=1 . 5 and a=30 
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Fig. 39 C T „ vs K for N=2 and ct=30 
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Fig. 40 C , vs K for N=2.5 and a=30 
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Fig. 41 C vs K for N=1.5 and a=30 
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Fig. 42 C _ vs K for N=2 and 



Single Cylinder 
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Fig. 43 C t t vs K for N=2.5 and 
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Fig. 44 C T . vs K for N=1.5 and a=30 
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Fig. 45 C . vs K for N=2 and 
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Fig. 46 C . vs K for N=2 . 5 and a=30 



Single Cylinder 





* 





* 



* 






2231 



* 






/ 



^ j 
* * 

/ 



/ 




o 

o 

\Q 

II 

3 




1-3 

u 



I t I I 1 I i 1 

co tn oj m -« m 3 * 

• ■ ■ 

ru — 



97 



Fig. 47 C T . vs K for N=1.5 and 



Single Cylinder 
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Fig. 48 C , vs K for N=2 and a=60 
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Fig. 49 C . vs K for N=2.5 and 
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Fig. 50 C, vs K for N=1.5 and a=60 
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Fig. 51 C vs K for N=2 and a=60 
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Fig. 52 C vs K for N=2 and 



Passive Cyl 
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Fig. 53 C vs K for N=1.5 and a=60 
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Fig. 55 C vs K for N-2.5 and a=60 
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Fig. 54 C T vs K for N=2 and a=60 
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Fig. 56 C . vs K for N=1 . 5 and a=60 
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Fig. 57 C tA vs K for N=2 and 
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Fig. 58 C T . vs K for N=2.5 and a=60 



